Normal metal to ferromagnetic superconductor tunneling 
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We study the point-contact tunneling between a normal metal and a ferromagnetic superconduc- 
tor. In the case of magnon-induced pairing the tunneling conductance is a continuous and smooth 
function of the applied voltage. For small values of the applied voltage the Ohm law holds. We show 
that one can obtain the magnetization and the superconducting order parameter from the tunneling 
conductance. In the case of paramagnon-induced superconductivity the tunneling does not depend 
on the magnetization. We argue that the tunneling experiment can unambiguously determine the 
correct pairing mechanism in the ferromagnetic superconductors. 
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The discovery of unconventional superconductivity 
caused an explosive growth of activities in various fields 
of condensed-matter research, stimulating studies of the 
basic mechanisms leading to this phenomenon. The most 
direct way to identify the Cooper pairs is from measure- 
ments of their spin susceptibility, which can be deter- 
mined by the Knight shift, from measurements of nuclear 
spin-lattice relaxation rate 1/Ti, probed by nuclear mag- 
netic resonance and nuclear quadrupole resonance, and 
by electron tunneling. 

In conventional superconductors, the quasi-particles 
form Cooper pairs in a spin-singlet state which has zero 
total spin. The existence of the gap in the quasi-particle 
spectrum leads to unusual properties of the systems: i) 
The specific heat decreases exponentially at low temper- 
ature, as opposed to the linear temperature dependence 
in the Fermi liquid theory. ii)The normal metal super- 
conductor tunneling experiments show that the electrons 
from the normal side of the junction can tunnel through 
and become an excited quasi-particle on the supercon- 
ducting side if the applied voltage is larger then the 
gap0. All these properties are well understood on the 
basis of the BCS theory of superconductivity. 

The discovery of superconductivity in a single crys- 
tal of UGe2^, URhGe^ and ZrZn2^ revived the in- 
terest in the coexistence of superconductivity and ferro- 
magnetism. The experiments indicate that the super- 
conductivity is confined to the ferromagnetic phase, the 
ferromagnetic order is stable within the superconducting 
phase (neutron scattering experiments), and the specific 
heat anomaly associated with the superconductivity in 
these materials appears to be absent. The specific heat 
depends on the temperature linearly at low temperature. 

At ambient pressure UGe2 is an itinerant ferromag- 
net below the Curie temperature Tc = 52K, with low- 
temperature ordered moment of /is — ^-'^I^b/U. With 
increasing pressure the system passes through two succes- 
sive quantum phase transition, from ferromagnetism to 
ferromagnetic superconductivity at P ^ 10 kbar, and at 
higher pressure Pc ^ 16 kbar to paramagnetism|20. At 
the pressure where the superconducting transition tem- 
perature is a maximum Tsc ~ 0.8/^, the ferromagnetic 
state is still stable with = "ilK. The survival of bulk 
ferromagnetism below Tac has been confirmed directly via 



elastic neutron scattering measurements |5|- The specific 
heat coefficient 7 = C jT increases steeply near 11 kbar 
and retains a large and nearly constant value |^. 

Specifically, U Ge2 has strong spin orbit-interaction 
that leads to an unusually large magneto-crystalline 
anisotropy with an easy magnetization axis along the 
shortest crystallographic axis. 

The ferromagnets ZrZn2 and URhGe are supercon- 
ducting at ambient pressure with superconducting crit- 
ical temperatures Tsc = 0.29K^ and T^c = 0.25K^ 
respectively. ZrZn2 is ferromagnetic below the Curie 
temperature Tc = 28. 5K with low-temperature ordered 
moment of /ij = 0.17/iB per formula unit, while for 
URhGe Tc = 9.5K and fi^ = 0A2fj,B- The low Curie 
temperatures and small ordered moments indicate that 
compounds are close to a ferromagnetic quantum critical 
point. 

We shall discuss two mechanisms of Cooper pairing 
in ferromagnetic metals: superconductivity induced by 
longitudinal spin fluctuations and magnon exchange 
mechanism of superconductivity'^. In the case of para- 
magnon induced superconductivity 7] the order param- 
eters are spin parallel components of the spin triplet. 
The theory predicts that spin up and spin down fermions 
form Cooper pairs, and hence the specific heat decreases 
exponentially at low temperature. The phenomenolog- 
ical theories 9] circumvent the problem assuming that 
only majority spin fermions form pairs, and hence the 
minority spin fermions contribute the asymptotic of the 
specific heat. The magnon exchange mechanism of super- 
conductivity was developed^ to explain in a natural way 
the fact that the superconductivity in UGe2 , ZrZn2 and 
URhGe is confined to the ferromagnetic phase. The order 
parameter is a spin anti-parallel component of a spin-1 
triplet with zero spin projection ("fl -|- H). The onset 
of superconductivity leads to the appearance of Fermi 
surfaces in the spin up and spin down momentum distri- 
bution functions. As a result, the specific heat depends 
on the temperature linearly, at low temperature. 

During the last years tunneling spectroscopy has been 
applied to identify the pairing symmetry. The most ba- 
sic idea of tunneling spectroscopy was first proposed by 
Bardeen|lOj who introduced the tunnel Hamiltonian ap- 
proximation for describing a tunnel junction. The con- 
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cept of the tunneling Hamiltonian|ll| became universally 
adopted for the discussion of tunneling in superconduc- 
tors. The idea is to write the Hamiltonian as a sum of 
three terms 



H — Hl + Hr + Hj 



(1) 



The first two terms Hj^ and Hjj are considered to be 
independent, expressed in terms of two independent sets 
of Fermi operators Ck,a, „ ™^ '^q,<yi ^9^0-1 where u = (t 
, 1). Tunneling is caused by the term Ht 



Ht 



^k.q'^q^.a-'^k.a 



(2) 



k,q,a 



The tunneling matrix element tk,q describes the transfer 
of the particle through the junction. In this letter we 
focus our attention on the point contact case, therefore 
it is assumed tk^q to depend only on the wave vectors on 
the two sides " fc" and "q" . The tunneling takes place over 
a very narrow span of energies near the Fermi surfaces 
of the "left" and "right" systems, that is why it is an 
adequate approximation to treat the transfer rate t^ q as 
a constant which is evaluated at kp and qp {tqp,kF = 

'^*qF,kp = 

We consider the tunneling process from a normal metal 
to a ferromagnetic superconductor. Hl is the Hamilto- 
nian of the system of free spin 1/2 fermions with disper- 

11, Hl = Y^^kCkaCka, and Hr is the 



sion efc 



2m 



k,a 



Hamiltonian for the ferromagnetic superconductor. For 
the magnon-induced superconductivity it has the form 



(3) 



^ Afe [dfetd-fei + dkid-k] + dtk^d+^ + rfl^^d^^ 

k 

while for paramagnon-induced superconductivity it is 



k,(T 



dfdk 



1 



dk^d-k^ 



- d+ 



(4) 

In Eqs.lEEIl EfeT = ^"M"-^ ^nd eki = ^-fJ.+H, where 
H is the energy of the Zeeman splitting, which is propor- 
tional to the low temperature ordered moment. There 
is a symmetry relation for the gap function Aj. which 
follows from the anti-commutation of spin 1/2 fermions 
A_fc = — Afc. The gap has the form A^ = Aq cos0fc|3,l3- 
The strong spin-orbit coupling in U Ge2 and U RhGe 
requires pseudo-spin techniaue|l2l| . The spin rotations, 
in that case, are accompanied by a rotation in momen- 
tum space. Hence if we consider tunneling process from 
a normal metal to a ferromagnetic superconductor, the 
tunneling Hamiltonian is invariant if and only if the tun- 
neling matrix element is a constant. This approximation 
is widely accepted, but in our case it is crucial for the 
applicability of the tunneling Hamiltonian. 




Next we calculate the current through the tunnel junc- 
tion. It is defined as 



(5) 



The tunneling current can be presented in the form 

I = etJ2 [< dg^^ 4a >< ^) - < Cka,d+„ >< (r, r)] 

k,q,a 

(6) 

where we have introduced the Green function in the 
Keldysh representation <C A,B~^ ^< (ri,T2) = 
i {B^ (T2) A{ti)) 13]. Let us consider first the magnon- 
induced superconductivity. We want to compute the tun- 
neling current to the lowest order (~ i^) in t[lj|- To that 
goal, we use the equations of motions of the Green func- 
tions in Eq. ® to cast the current into the form 



/ = et 



where 



k,q,cr 



/duj 



{q,Uj) [G^r {k,Lu)-G,a (fc, C^)]] 



(7) 



Et^(g,a;) = ulAi^iq,^)+v^qA+^iq,ij), (8) 
= vlA'l^{q,uj) + ulA2f,{q,uj), fi = r,a,< . 

Here 



ui = - I 



2 1 2 



(9) 



The expressions for the retarded/advanced Green func- 
tions are given as follows 

G,r/a{k,u;) = (w-efe±zO+)"' (10) 
Ai,/a{q,uj) - {Lj-Eiq±i0+y\l^l,2 

and cj) = -Ai^^^{q,-iL;). The quasiparticle en- 

ergies are 



Eiq = -H- ,/e2 + A2, 



H 



A2. 



(11) 
(12) 



^2q y . — 

The distribution Green function is defined as 

G,< (fc, LO) = -/fd(^) [Gar (fc, ^) - G„a [k, u)] , (13) 

/fd('^) is the Fermi-Dirac function and likewise for 
Ai^,l — 1,2. Next we convert the sums over the wave 
vectors into integrals over the corresponding energies by 
introducing the density of states for the left-hand and 
the right-hand side of the junction. We consider the sim- 
plest possibility of constant density of states, that is, we 
do not take into account the finite band-width effects. 



This choice aUows us to obtain analytical results in the 
zero-temperature case. The expression for the tunneling 
current is obtained in the form: 
I) If Ao < 2H 
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II) If Ao > 2H 
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where /o = ■net^ p-^'' p-^'' /^q,x± 



d) , < < Ao - 

d)-/(a:+), /^o-H<eV<H 
d)-/(x+), i7<ey<Ao + 7J 
-/(.T^+)-/(a;-),el/> Ao + i? 

(14) 

0<ey<i/ 

•^), H<eV<Ao-H 
■l) - f {x+) ,Ao-H <eV <Ao + H 
■l)-f{x+)-f{x^), eV>Ao + H 

(15) 

{eV±H)/AoJ{x) = 
is the density 



x^ tan~^ ^/x'^ — 1 — s/x^ — \ and p)^' 
of states in the left-hand/right-hand side of the junc- 
tion. Also, we calculated the dynamical conductance 
g = dl / dV. The expression for it is given by 
I) If Ao < 2H 
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0<eV <Ao-H 
Ao~ H <eV < H 
H <eV <Ao + H 
gix-),eV> Ao + H 



(16) 



II) If Ao > 2H 



9_ 

90 



Ao ' 
TreV 
Ao ' 



Ao ' ^ ^ , ^ 

^-5(a;+),Ao-i/<ey < Ao 
l!^^g{x+)~g{x-), eV>A, 



0<eV <H 
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+ H 
+ H 



(17) 



where go — 2iTe'^t^ p"^^ p^^'^ and g {x) = a; tan ^ \/ x'^ 



1. 



The results for both the tunneling current and the dif- 
ferential conductance are shown in Fig. 1 and Fig. 2. 
The most important feature of the tunneling conductance 
is that it is a continuous and smooth function of the ap- 
plied voltage and does not show any discontinuities at 
the gap edges that are characteristics of the tunneling in 
which conventional superconductors are involved. Also, 
the tunneling current is non-zero when eV < Aq. One 
observes that for small values of the applied voltage / 
depends linearly on it, that is, the Ohm law holds. This 
behavior is attributed to the existence of a Fermi sur- 
face in this model'^. Another important feature to be 
pointed out is that the magnetization which is included 
in the Zeeman splitting energy H can be measured in 
a tunneling experiment. When Aq < 2H, the differ- 
ential conductance is constant for eV < Aq — H, the 
minimum is located at eV = H and the maximum is at 
eV = Aq + H [Fig. 2, the solid curve]. In the case of 
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FIG. 2: The differential conductance as a function of the ap- 
plied voltage in the case: i) magnon-induced superconduc- 
tivity for Ao = l.^H (solid line) and Ao = 3-ff (short-dashed 
line) ; ii) paramagnon- induced superconductivity (long-dashed 
line) . 



Aq > 2H, g exhibits ohmic behavior for eV < H and the 
two maximums are at Aq — H and Ao -I- H [Fig. 2, the 
short-dashed curve]. Thus, for any value of the micro- 
scopic parameters one can deduce the values of both Ao 
and H from the tunneling conductance. 

Let us consider now the paramagnon-induced super- 
conductivity. In this case, the tunneling current is a sum 
of two terms / = In+Isc where In is the normal current 
due to the unpaired band of down-spin electrons and Isc 
is the current of the superconducting up-spin electrons. 



The explicit expressions for them are written as 



In = ei^^ / ^[[D^riq.io)-D^a{q,Lo)]G^<ik,io) 



k,q 



(IS 



and 



duj , 



Isc = et^J2 I ^[P.(g,^)-Sa(g,a;)]G|<(fc,c^) 



(19) 



where 



^ u'^gA^{q,uj) + VgA+{q,uj), fi^r,a,<. (20) 

The retarded/ advanced Green function are given 
by Dir/aiq,^^) = (t^ - e<;i ± «0+) \ A^/^ (g, w) = 
(u — Eq ± iO'^) and the distribution Green functions 
are defined as above. The quasiparticle energy is 



E„ 



(21) 



Also, 



"^2 1 E„ 



^ 2\ E, 



(22) 



With the same assumptions as above the tunneling cur- 
rents are obtained in the form 

In = 2/oey/Ao (23) 



Isc = 2/o 



Unlike the case of magnon-induced superconductivity the 
current does not depend on the Zeeman splitting energy. 
The corresponding expression for the differential conduc- 
tance is cast into the form 



o 
cn 




go 



(25) 



The results for the tunneling current and the differen- 
tial conductance are shown in Fig. 1 and Fig. 2 (the 
long-dashed- line), respectively. One sees that g has no 
ohmic behavior and has only one maximum located at 
Aq. Thus, a tunneling experiment can very easily deter- 
mine the correct pairing mechanism in the ferromagnetic 
superconductors. 

Point-contact tunneling between normal metal and su- 
perconductor is one of the best probes for analyzing the 
energy gap of superconductors. Measurements of tun- 
neling current yield valuable information on the symme- 
try of the order parameter which in turn is essential for 



eV/Ao 

FIG. 4: Bogoliubov quasiparticle contribution to the differ- 
ential conductance in the case of magnon-induced supercon- 
ductivity for Ao = 3H: quasiparticle with energy Eiq [Eq. 
(llH l (short-dashed line) and quasiparticle with energy E2q 
[Eq. GU] (solid line). 



understanding the mechanism of superconductivity. Our 
results differ from what we expect for the conventional su- 
perconductors in many aspects. Perhaps the most strik- 
ing difference is the nonzero conductance inside the su- 
perconducting gap. The tunneling conductance is a con- 
tinuous and smooth function of the applied voltage. For 
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small values of the applied voltage / depends linearly on 
it, that is, the Ohm law holds (Figl). The main focus 
in the present paper is the evidence that the origin of 
this behavior is the existence of a Fermi surface in the 
superconducting phase. It was experimentally observed 
0, 01 that the specific heat depends on the temperature 
linearly at low temperature. This means that the super- 
conducting state is strongly gapless with a Fermi surface. 
Our theoretical prediction suggests that the tunneling ex- 
periment could give another independent verification of 
the surviving of the Fermi surface in the superconducting 
state. 

We consider two pairing mechanisms - magnon - and 
paramagnon- induced superconductivity 0, the 
case of magnon induced superconductivity the differen- 
tial conductance exhibits ohmic behaviour at low voltages 
for any value of the microscopic parameters, and has two 
extreme points which determine the gap and the Zee- 
man splitting. Unlike this, the tunneling current, in the 
case of paramagnon-induced superconductivity, does not 
depend on the Zeeman splitting energy, the differential 
conductance has no ohmic behavior, and the only local 
maximum determines the gap. 

The contribution of the Bogoliubov quasiparticles to 
the differential conductance is shown in Fig. 3 (for 
Aq = 1.5H) and in Fig. 4 (for Aq = 3H) in the case of 
magnon-induced superconductivity. From Eqs. H11I12|) 
it can easily be established that there is a gap Ao + H in 
the quasiparticle energy Eiq while the quasiparticle with 
energy E2q is a gapless excitation with a Fermi surface 
As a result, the contribution of the first quasiparticle to 
the differential conductance shows typical behaviour well 
known for gapped superconductors [Figs. 3 and 4, the 
short-dashed line] and all important features in g we have 
already pointed out are due to the second quasiparticle 
[Figs 3 and 4, the solid lines]. Namely, the ohmic part in 
the differential conductance at low voltage bias reflects 
the gapless nature with a Fermi surface of this quasipar- 
ticle. The ferromagetic superconductors are anisotropic. 
The different symmetries of the order parameter corre- 
spond to different anisotropics, respectively to different 
Fermi surf aces. As a consequence, calculating the tunnel- 
ing current and the differential conductance and averag- 
ing over the Fermi surface we obtain different expressions 
for the different Cooper pairing mechanisms. 

It is important to compare our results with the results 
in the case of isotropic s-wave superconductors. At zero 
magnetic field the energies of the Bogoliubov quasipar- 
ticles are degenerate and there is only one peak in the 
quasiparticle density of states. Correspondingly, there is 
only one peak in the differential conductance. The ef- 
fect of the magnetic field is to lift this degeneracy. The 
quasiparticle peak in the density of states splits in two 
peaks and as a result two peaks appear in the differential 
conductance with each peak comin g fr om a quasiparticle 
with a given spin (up or down) [T5Lll6j] . 



In the case of magnon-induced superconductivity the 
quasiparticle spectrum is not degenerate too. The two ex- 
trema in the differential conductance result from a com- 
plicated averaging over the Fermi surface and reflect the 
symmetry of the superconducting order parameter, re- 
spectively the mechanism of the superconductivity. 

In the case of paramagnon-induced superconductivity 
there is one gapped excitation H21|l which gives the only 
peak in the differential conductance. The other excita- 
tion is a free spin-down electron which is responsible for 
the linear in the bias voltage part of the tunneling current 
(El- 
It is known that in the case of anisotropic pairing zero- 
energy bound states can exist [l7j . They will modify 
the low-voltage behaviour of the differential conductance. 
The inclusion of these states would require modification 
of our approach. The results we have obtained will be 
changed at low voltages. However, the non-zero differen- 
tial conductance at zero voltage is a direct consequence of 
the existence of a Fermi surface in the ferromagnetic su- 
perconductors. This behaviour will not be changed quali- 
tatively when the zero-energy bound states are taken into 
account. The rest of out findings should apply in this 
case because they refer to features in the differential con- 
ductance at voltages of the order of the superconducting 
order parameter. The evident difference between the tun- 
neling processes involving magnon-induced superconduc- 
tivity and those involving paramagnon-induced super- 
conductivity suggests that tunneling experiment could 
be crucial for understanding the mechanism of ferromag- 
netic superconductivity. 

The ferromagnetic superconductivity was seen only in 
samples with small normal-state residual electrical resis- 
tivity pq — 2/iricm 3, 4]. This means that the impurity 
concentration is very low. Therefore, if the impurities 
are taken into account they will only slightly modify the 
low-voltage behaviour of the tunneling differential con- 
ductance. 

The tunneling differential conductance in the case of 
s-wave superconductors in an external magnetic field has 
been measured long ago It shows different behav- 

ior for different values of the magnetic field [Fig. 1 in 
Ref.|l5j]]. that is, the tunneling current measurements 
are sensitive to the quasiparticle spectrum of the sys- 
tem in the magnetic field. This fact demonstrates that 
the experiments can in principle distinguish the different 
behavior we predict for the two pairing mechanisms. 

Our results differ, as well, from the calculations of the 
tunneling conductance in the case when spin-trjplet su- 
perconductivity of different kind is considered 18]. This 
demonstrates that by means of tunneling experiments the 
symmetry of different p-wave superconductors can be es- 
tablished. This is an additional support for our efforts to 
distinguish the magnon and paramagnon mechanisms of 
FM superconductivity. 
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